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, Abstract. Magnetic stresses cohimate protostellar winds into a com- 

mon distribution of force with angle. Sweeping into the ambient medium, 
such winds drive bipolar molecular outflows whose properties are insen- 
sitive to the distribution of ambient gas and to the details of how the 
. wind is launched, and how its intensity varies over time. Moreover, these 

properties are in accord with the commonly observed features of outflows. 

This model is simple enough to permit a quantitative study of the 
feedback effects from low-mass star formation. It predicts the rate at 
' which star-forming gas is ejected by winds, and hence the efficiency with 

\ which stars form. Applied to individual star formation, it relates the stel- 

lar initial mass function to the distribution of pre-stellar cores. Applied 
\ to cluster formation, it indicates whether the resulting stellar system will 

remain gravitationally bound. 
, Using the energy injection and mass ejection implied by this model, 

' we investigate the dynamical evolution of a molecular clump as a stellar 

Q^. cluster forms within it. This depends critically on the rate at which tur- 

bulence decays: it may involve equilibrium star formation (slow decay), 
overstable oscillations, or collapse (fast decay). 
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■ Stars, even low-mass stars, are not born quietly. In the process of accreting gas 

through an accretion disk, each star returns a fraction of the inspiralling material 
to the interstellar medium at speeds comparable to the escape velocity from the 
stellar surface. These powerful winds, observable as jets of ionized gas, sweep 
the ambient material into opposed streams known as bipolar molecular outflows. 
Because of their intensity, protostellar winds pose severe challenges to those who 
wish to model star forming regions. On a small scale, the wind from a forming 
star is capable of stripping gas from its parent molecular core even as that core 
feeds material to the star. On a larger scale, the outflow driven by a star's 
wind can punch through a massive molecular clump that is forming a stellar 
cluster. There is now plentiful observational evidence that these interactions 
are common, even ubiquitous, within star-forming clouds. Moreover, there is 
evidence to support the argument (Norman & Silk 1980) that the feedback from 
protostellar winds is the primary driver of turbulence within molecular clouds, 
and may be responsible for establishing Larson's (1981) line width-size relation 
in these clouds (McKee 1989). 
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Because these phenomena have profound consequences for the formation of 
stars and star clusters and for the evolution of molecular clouds, it is important 
to quantify the ways in which protostellar winds impact their molecular envi- 
ronment. In the past, models of feedback in star formation have either ignored 
the effects of low-mass stars (e.g., Elmegreen 1983) or considered only spheri- 
cally symmetric protostellar winds (e.g., Nakano, Hascgawa, & Norman 1995), 
despite the fact that these winds are clearly jet-like. More sophisticated models 
of protostellar outflows (such as those of Shu ct al. 1991; Masson & Chcrnin 
1993; Li & Shu 1996) have not been applied to the basic questions of feedback. 
Partly, this has been caused by a lack of consensus (sec Richer ct al. 2000) as to 
which of the several classes of models for outflows best describes their dynamics. 

We shall argue that there is a single model for protostellar outflows that 
is both motivated by theories of protostellar winds and validated by its ability 
to explain some of the common, otherwise mysterious features of protostellar 
outflows. Although this model cannot explain all the details of these flows, 
it provides a stable platform from which to launch an exploration wind-cloud 
interactions. We shall then use this model to address basic questions of star 
formation: the stellar initial mass function (IMF); the efficiency of star formation 
in clusters; and the dynamical evolution of star-forming gas during stellar cluster 
formation. 

The work presented here has been conducted in collaboration with Chris 
McKee and appears or is intended to appear in: Matzner (1999), Matzner & 
McKee (1999a), Matzner & McKee (2000), and Matzner, Bertoldi, & McKee 
(2000). 

2. Common Properties of Bipolar Molecular Outflows 

Protostellar molecular outflows typically share several common features, as dis- 
cussed by Masson Sz Chernin (1993) and Lada Sz Fich (1996). These are: a 
roughly linear position- velocity (PV) diagram (a "Hubble law"); a lack of reced- 
ing material in the approaching lobe, and vice versa; and a power law distribution 
of mass with velocity: dM /vohs oc v^^^, where fobs is the line-of-sight velocity 
relative to the systemic velocity and T is typically about -1.8. Of these, the 
power law distribution of mass with velocity is the most difficult to explain, as 
it often holds quite well over a factor of five or ten in fobs (e.g., Lada & Fich 
1996). 

In order to construct a theory for the feedback effects of protostellar winds 
it is necessary to develop a model for these outflows, which represent the inter- 
action between winds and the ambient gas. To be tenable, such a model must 
confront the common properties of outflows listed above. Moreover, it must 
also be compatible with the ways in which protostellar winds are launched and 
collimated, a topic we now address. 

3. Structure and Intensity of Protostellar Winds 

Protostellar winds are launched centrifugally from the inner regions of proto- 
stellar accretion disks (Blandford & Payne 1982). This is made possible by the 
presence of intense poloidal magnetic fields, which act as rigid tubes anchored in 
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the disk, along which gas is flung away. Once the wind has traveled sufficiently 
far, the field weakens and gas begins to travel ballistically rather than being 
forced to corotate with the disk. Since it is still stuck to the same field line as 
when it was launched, and since that field line stretches back to the rotating 
disk, the magnetic field is wound into a tight spiral. By the Biot-Savart law, the 
tightly- wound field implies that the wind encloses a current. Magnetic stresses 
associated with current within the wind can generate forces that push the wind 
toward or away from its axis. These forces arc important for the collimation of 
the wind on scales ( ^ 1 AU) comparable to the launching region of the accre- 
tion disk, but on the larger scales of molecular outflows (~ 0.1 pc), the wind 
relaxes into transverse pressure balance. Therefore, the wind's magnetic flux 
lines tend toward a force-free state in which there is no current within the wind: 
the enclosed current lies entirely along the axis, rather than threading the wind 
itself. The field thus settles into the force-free state i?^ = 21{r)/ (cw), where 
w = r sin 9 is the cylindrical radius and the current /(r) is a very slowly varying 
function of r. (There must be currents, and hence forces, in the wind for I{r) 
to vary at all: these effect a continuing gradual collimation of the wind.) 

This force-free distribution of toroidal magnetic field is sufficient to deter- 
mine the ram pressure of the wind at large distances from its source. For, the 
magnetic field lines must corotate with the disk in a steady state. This, along 
with the fact that the wind will expand to fill the entire solid angle available 
to it, allows one to trace the field lines from their origin on the disk to their 
destination on the sky (Ostrikcr 1997; Matzncr 1999). Then, the conservation 
of mass and energy along each streamline determines the ram pressure of the 
wind in each direction. In the case of a wind whose density just above the disk 
varies with disk radius zuq as pu,Q oc vjq , and whose field lines expand signifi- 
cantly (tu ^ zuq) according to the force-free distribution of the field, the wind 

ram pressure varies as Pwv'^ oc C{r)'nj^ '^^'^w~'^ (Ostriker 1997). Now, the disk 
radius tuo is generally restricted within a very narrow range compared to the 
wind radius w. for instance, in the model of Shu et al. (1995), vjq takes only 
a single value whereas vj extends from a few AU to ~ 0.1 pc. For this reason, 
the ram pressure of the wind is very close to Pw^w w"^ in practice (Matzner 
& McKee 1999a). This force distribution, which is characteristic of the x-wind 
model (Shu et al. 1995), is thus common among hydromagnetic winds. 

In reality, variability and instabilities are likely to broaden the wind force 
inside some angle ^o; therefore, we consider the smoothed distribution 

^ m^Vy, 1 1 

P^""^ - 47rr2 ln(2/0o) (sin Of + (sin ^o)^ ' ^ ' 

In the model for outflows presented below, observations imply 9q ~ 0.01. 

The fact that hydromagnetic winds tend to PwVw oc (sin^)~^ is consis- 
tent with the theoretical expectation (Matzner & McKee 1999a) that the fast 
Alfvenic Mach number, B'^ / (An p^v'^) , is roughly constant. This expectation 
follows from the fact that this ratio is unity near the source, and varies loga- 
rithmically with distance along each streamline. The (sin^)~^ force distribution 
implies that the cumulative wind momentum within some angle, Pw{< &) = 
2TTr'^p^Vw sin 6' dd'dt, varies as ln(^/0o) for a broad range of angles, and this 
property will be useful for understanding outflows and their effects. 
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Figure 1. Mass- velocity relations for L1551 {triangles, Moriarty-Schieven 
& Snell 1988), NGC2071 {squares: Moriarty-Schieven, Hughes, & Snell 1989), 
and NGC2264G {diamonds: Lada & Fich 1996), fit with outflow model solid 
lines). The inclination has been chosen to make the spectral break near the 
lowest observed velocities: 50° for L1551, 40° for NGC2071 and NGC2264G. 
Filled symbols: blue lobes; open symbols: red lobes. Velocities are relative to 
vo, a center velocity, chosen to maximize the symmetry between each pair of 
lobes. Curves have been offset vertically for clarity. 

A protostar's wind lasts as long as it accretes material through its disk, 
which is essentially the free-fall time of the unstable core from which it formed, 
typically ~ 10^ yr. If the wind's mass is a fraction of the star's mass, 
then the wind momentum is fwVw'^-k- We estimate that fwVw — 40 km s""^, 
roughly consistent with the observational estimate of Richer et al. (2000); for a 
wind velocity ~ 200 km s^^, this implies fw ~ 1/5, between the theoretical 
predictions of Shu et al. (1988) and Pelletier and Pudritz (1992). We expect 
fvjVw to be approximately independent of stellar mass, because Vyj scales roughly 
with the stellar escape velocity, which itself is regulated by deuterium burning 
during accretion (Stabler 1988). 

4. Bipolar Molecular Outflows from Hydromagnetic Protostellar Winds 

The protostellar winds described above are essentially radial, both because they 
expand to fill the available angle, and because the wind coasts after achieving 
force balance. At wind speeds typical of low-mass protostars, the shocks that 
separate the wind from the surrounding gas are radiative (Koo &: McKee 1992), 
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so the swcpt-up shell conserves momentum in each direction. For a steady 
wind and an ambient medium whose density varies as Pa oc Q{9)r~^, Shu et al. 
(1991) showed that these ingredients reproduce the Hubble law for outflows. 
However, Matzner &; McKee (1999a) demonstrated that neither a steady wind 
nor pa oc is necessary for this conclusion: as long as the ambient medium is 
relatively featureless (e.g., a power law) on the scale of the outflow, the outflow 
motion is self-similar even if the wind intensity varies. For radial motion this 
implies the Hubble law, since v{9,t) oc r{9,t)/t. The collimated wind force 
distribution Pw^w oc (sinO)^^ leads to outflows that are highly elongated and 
have bowshock-shaped tips of width r^o- This morphology is thus consistent 
with a "wide-angle" driving wind. A collimated outflow with v oc r lacks blue 
material in its red lobe and vice versa, consistent with observations. 

That leaves only the power-law distribution of mass with velocity, dm/dvobs 
"^olas^' explained. Masson & Chernin (1992) argued that the radial outflow 

model of Shu et al. (1991) was inconsistent with this relation, but they only 
considered angular distributions Pwv"^ oc (cos 6')"''' for some /3. Their results 
do not apply to hydromagnetic winds, for which oc (sin^)~^. Indeed, Li 

& Shu (1996) demonstrated that the outflow driven by a steady x-wind into a 
magnetically-flattened core, with pa oc Q{6)r~'^, could produce dm/dvo\ys oc v^^^. 
I shall now show that this conclusion is much more general (see also Matzner & 
McKee 1999a), and is essentially independent of the ambient medium. 

The outflow shell expands at the rate dictated by momentum conservation; 
since pwV^ on for small 9, and since one expects the ambient medium also to 
be a power law of 9, the rate of expansion obeys a power-law relation v{9) oc 9~^ 
quite generally. The power x can be computed, but it is not important for the 
present argument. Now, hydromagnetic winds satisfy the cumulative momentum 
distribution Pw{< 9) oc ln(0/0o) for a wide range of angles, as described in §3. 
But since v oc 9~^, momentum conservation [psheii(< = Pw{< ^)] implies a cu- 
mulative momentum distribution with outflow velocity Psheii(> u) oc ln('y/vmax), 
where fmax is the rate of expansion along the wind axis. This, in turn, implies 
that dm/dv = (i^Pshcii(> v)/dv'^ oc and, since i;obs oc v for an elongated 
outflow, dm/dvohs oc v^^^, i.e., F = —2. Notice that essentially nothing about 
the ambient medium entered into this argument; therefore, F ~ — 2 is generic. 
Figure 1 demonstrates that models with this scaling can reproduce the mass- 
velocity distributions of real outflows. 

The commonly-observed features of molecular outflows are thus the natural 
products of momentum-conserving shells driven by hydromagnetic protostellar 
winds: they are insensitive to the details of the wind's launching region, of vari- 
ations in the wind's intensity, and of the ambient density distribution. Although 
this model remains to be tested in detail, it matches observation well enough to 
deserve attention as a mechanism for feedback. 

5. The efficiency of star and star cluster formation and the IMF 

The intensity of a protostellar wind along its axis ensures that within some angle, 
gas will be blown away. Because nearby gas is typically accreting onto the star 
itself or in the process of forming other stars, this process limits the efficiency 
with which individual stars and multiple stellar systems can form (see Matzner 
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& McKee 2000, for a more detailed discussion). Gas dispersal by protostellar 
winds is more ubiquitous and less violent than disruption by massive stars. 

To estimate the amount of gas (Mej) lost per outflow, first note that outflows 
propagate by conserving momentum in each direction. To be ejected, material 
must travel faster than the escape velocity Vesc of the system. The momentum of 
the escaping gas is at least M^jVesc] therefore, Mejt^esc < fwVwfn*- In reality, mo- 
mentum is lost in those directions where the wind is too weak to eject anything, 
and wasted in those where the ejected gas travels faster than v^sc- Since the 
wind momentum is distributed logarithmically with angle, the fraction available 
to eject material at around v^sc is roughly ln(l/0o); in fact, ln(2/^o) turns out to 
be a better estimate. Moreover, an outflow's momentum is reduced by a factor 
Cg (typically of order unity) due to the action of gravity as it traverses the cloud. 
Defining the efficiency parameter X and the efficiency per star e, 

X = Cg ln(2/6'o)ucsc/ ifwVw), £ = rri^l (m^ + Afcj), (2) 

Matzner & McKee (2000) find that for the formation of an individual star from 
a collapsing molecular core, 

e-' (2X)-i + [(2X)-2 + (1 + /^)2]i/2. (3) 

and for the formation of a star within a larger cluster-forming clump, 

e"^~l + (2X)-^ (4) 

In both cases, £ oc X for X <C 1 (strong winds) and e — > 1 for X ^ 1 
(weak winds). The core and clump cases differ for a number of reasons. In the 
formation of a single star, unlike in multiple star formation, we can safely assume 
that all of the core mass is either accreted or ejected. Clumps have higher values 
of 

Vesc (which raises X), but their density profiles are shallower (which lowers 
X through Cg); these eff'ects roughly cancel. Finally, mass is ejected at angles 
much further from the wind axis in the core case; for this reason, the ejected 
wind mass is included in (|3|) but neglected in 

For the formation of an individual star, (^) assumes that the collapsing core 
is spherically symmetric. However, protostellar cores are likely to be partially 
supported by magnetic fields, which will cause them to flatten along the magnetic 
axis. Since the protostellar wind is likely to share this axis, magnetic flattening 
segregates material away from the wind axis, reducing the amount that will 
be ejected and increasing the efficiency. Employing the models of flattened, 
magnetized cores presented by Li & Shu (1996), Matzner & McKee (2000) find 
that equation (|^) remains valid if X is replaced by (1 + Hq)^/'^X, where 1 < 
(1 + Hq) ^ 2 is the factor by which magnetic support increases the core's mean 
density. In the limit of a completely disk-like core {Hq ^1), only the wind mass 
would be lost [e — >■ 1/(1 + fw)]- Individual stars form at about 30% efficiency 
if cores are spherical {Hq = 0), or ~ 75% if they are significantly fiattened 
{Hq = 1), as shown in the left panel of figure 2. 

For the formation of a star within a stellar cluster- forming clump, the 
anisotropy of the region is not likely to correlate with the wind axis. But, 
(^ assumes that each star forms at the center of a spherical distribution. Com- 
paring again simulations in which stars form in a distributed manner throughout 
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Figure 2. Left panel: Core disruption and the efficiency of single star for- 
mation. The efficiency e of protostellar core collapse (equation p|), versus the 
efhciency factor X. Curves are plotted for spherical {Hq = 0), moderately 
flattened {Hq = 0.4) and significantly flattened {Hq — 1) cores. Circles indi- 
cate the values of X for thermally supported (T = 10 K) cores; turbulence 
can only increase X above these values. Right panel: Estimates of the in- 
stantaneous star formation efficiency s in stellar cluster formation. Dashed 
lines: analytical estimate, equation (Q). Solid lines: Numerical evaluation of 
the same quantity in a protostellar clump with the structure p (x and in 
which stars form at the local ambipolar diffusion rate. Dotted lines indicate 
the mean extinction of the clump. Above the dash-dot line, stars more mas- 
sive than 20 Mq are likely to form; our estimates are therefore only rough 
upper limits in this region. The analytical expression relies on assumptions 
that break down outside of the upper and lower thin lines. 



a cloud, Matzner & McKee (2000) find that (^ is an excellent approximation 
nevertheless, as shown in the right panel of figure 2. Equation (^) predicts 
efficiencies 30% ^ £ ^ 50% for the conditions typical of low-mass stellar cluster- 
forming regions, like those in Orion B (Lada 1992): each star expels one to two 
times its own mass from the clump. It is difficult to compare E with its obser- 
vational analogue, the ratio of stellar to total mass, but the two quantities are 
similar in a few well-studied cluster- forming regions (Matzner &: McKee 2000). 

Using the efficiency of individual star formation, we can estimate the de- 
gree to which the stellar IMF differs from the mass function of pre-stellar cores 
(Nakano et al. 1995), a topic that has recently come under observational scrutiny 
(Motte, Andre, &: Neri 1998). This comparison, discussed in detail by Matzner 
& McKee (2000), involves a consideration of how much X, and thus e, differs 
among cores of different masses. We find that the variation of e is quite subtle, 
and that the slope of the IMF is thus only slightly flatter than the slope of the 
core mass function. 

It should be noted that equations (P) and (Q) predict efficiencies significantly 
higher than those found by Nakano et al. (1995), and that this is a direct result 
of the collimation of protostellar winds. By this process (in contrast to the effect 
of a massive star), mass is lost continuously as stars form. The efficiencies found 
above are consistent with the formation of bound clusters (Mathieu 1983), which 
raises the question of why so few actually do remain bound (Lada &: Lada 1991). 
This could be due to a violent event associated with massive star formation; 
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alternatively, clusters could be unbound by motions of their parent clumps, a 
topic we address below. 

6. Dynamical Evolution of Clumps During Cluster Formation 

Embedded stellar clusters are found exclusively within the most massive clumps 
inside molecular clouds. Since these clumps are massive enough to be confined 
by their own self-gravity despite the weight of the molecular cloud above them 
(McKee 1999), we can gain insight into the production of a stellar cluster by 
approximating a clump as a distinct reservoir of gas, and examining its dynam- 
ical evolution while a stellar cluster forms within it. This is the subject of an 
ongoing investigation (Matzner 1999; Matzner & McKee 1999b; Matzner et al. 
2000), and the results reported here should be considered preliminary. 

Like any self-gravitating gaseous system, a cluster-forming clump must sat- 
isfy the virial theorem for self-gravitating gas (here in Eulerian form; McKee &: 
Zweibel 1992): 

^I = 2{T-To)+W + M-~ |^(pvr2) • dS, (5) 

where / ~ MciR'^i/2 is the trace of the clumps's moment of inertia tensor; W 
and A4 are its gravitational and magnetic energies; the kinetic energy of its 
internal motions is T; Tq is an energy associated with its confining pressure; 
and the last term accounts for gas that may leave the system. Equation (|5|) has 
been used to identify equilibrium states of star-forming molecular clouds (McKee 
1989; Bertoldi & McKee 1996), but it can also be taken as a dynamical equation 
if not taken to equal zero. In order to accomplish this, each of the energetic 
terms must be approximated as a function of four dynamical variables: the 
clump's mass Md and the mass of its embedded cluster, its radius Rd, and 
its internal velocity dispersion Vrms- To simplify matters, the results presented 
here assume a spherical, magnetically supercritical {M^i = 2M$) clump with a 
p oc density profile and no rotation; any stars are assumed to move along 
with the gas. 

To complete our description of the system, it is necessary to specify the 
evolution of Md, M^, and Urms- All three are affected by star formation and 
the effects of protostellar outflows on the gas. For this, we use the theory of 
McKee (1989): the local rate of gas conversion into stars (tg^) is identical to 

the local ambipolar-diffusion rate (^ad)' because the formation of protostellar 
cores requires the assembly of magnetically supercritical condensations. Because 
of ionization by external FUV photons, star formation is inhibited within four 
visual magnitudes of the surface; interior to that, it proceeds at about a twelfth 
of the local free-fall rate. As stars form, their mass is added to and subtracted 
from Mci; also subtracted is the mass of gas that is ejected according to equation 

For a massive molecular clump, t^rms is dominated by non-thermal turbu- 
lent motions. The turbulent energy decays at a characteristic rate that can be 
parameterized by (dw^rms/'^i) decay = ~'>^Tms /('its), where tg is the free-fall time 
of the clump. Therefore, turbulence must be replenished every r] free-fall times; 
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Figure 3. The evolution of a 500 Mq molecular clump as a stellar cluster 
forms within it. All simulations start from a state of equilibrium star forma- 
tion; they differ in the rate of turbulent decay {rits)~^, and in whether star 
formation is taken to be continuous {solid lines) or individual {dashed lines). 



we expect 1 r/ < 10, where numerical experiments (e.g., Stone, Ostriker, & 
Gammie 1998) favor r] ^ 1. Turbulent energy can derive from two sources. For 
one, it can be exchanged for other energies of the system, through work done 
during compressions or expansions of the cloud; for this, we use the analytical 
results of Zweibcl &: McKee (1995) and McKee & Zweibcl (1995) to estimate 
the energetics of turbulent, magnetized gas. This process does not, however, 
prevent the cloud from collapsing. Alternatively, turbulence may be stirred up 
by protostellar outflows. Since these conserve momentum, it is most convenient 
to consider the clump's turbulent momentum (Mcifrms) enhanced an amount 
4'wPw by a wind of momentum p^j-, here, 4>w ~ 0.6 is the fraction of the wind 
momentum that does not escape the clump (Matzner 1999). 

Putting these elements together, it is useful to consider two approximations: 
one in which energy is added continuously according to the mean star formation 
rate within the clump, and another where stars form randomly at the current 
star formation rate, are chosen individually from the IMF, and affect the cloud 
instantaneously. These approaches differ for clumps of about five hundred solar 
masses or less, in which only ten to twenty stars form per tg. In the continuous 
approximation it is possible to identify equilibria (Norman &; Silk 1980; McKee 
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1989) in which turbulent decay is just offset by injection from outflows. These 
are stable, in the sense that a compressed cloud's turbulence is enhanced faster 
than it decays. But they can also be overstable, as energy in an oscillation can 
grow during compressive phases. 

Perturbations due to individual forming stars will stimulate oscillations of 
the system, whose response depends critically on the rate of turbulent decay. 
Figure 3 shows the effect of changing r] for a clump of SOOM©: when turbulence 
decays slowly (r/ = 10), a cluster- forming cloud will remain very close to its 
equilibrium state, whether or not star formation is taken to be continuous. For 
moderate tiubulent decay (// = 3), individual star formation seeds overstable 
oscillations about equilibrium: stars form in bursts. For fast turbulent decay 
{t] = 1), a continuous simulation can remain near its equilibrium, but one with 
individual star formation cannot consistently replenish the lost turbulent energy 
every dynamical time; the system collapses. Although the approximations used 
in these simulations (e.g., spherical symmetry and no rotation) affect these out- 
comes, it appears that the duration and constancy of star formation are sensitive 
to the rate of turbulent decay. 



7. Conclusions 

Observations of star-forming regions and of individual accreting protostars have 
made it increasingly clear that stars interact violently with their parent clouds as 
they form. To model this phenomenon requires an understanding of protostellar 
winds, and of the bipolar molecular outflows they drive into the surrounding 
gas. We have seen that the typically observed features of molecular outflows 
can be explained as the natural product of magnetically-collimated winds and 
radiative shocks. This yields a model for outflows which can be used to estimate 
how much mass is ejected, either from an individual collapsing core or from a 
massive clump in the process of forming a stellar cluster. The estimates of star 
formation efficiency (30% — 50%) presented in §5 are much higher than those of 
Nakano et al. (1995), because they account for the collimation of protostellar 
winds. 

Lastly, the feedback effects of protostellar outflows can be incorporated 
into simple dynamical models for the process of stellar cluster formation within 
massive molecular clumps. The investigations described in §6 show that the 
dynamics of stellar cluster formation are sensitive to the rate at which turbulence 
decays, and this will have observational consequences. 

In the future, large-scale numerical simulations of star formation must in- 
clude protostellar outflows, for without them such simulations are incomplete. 
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